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Abstract 

Structural and kinetic aspects of 2-D irreversible metal deposition 
under potentiostatic conditions are analyzed by means of dynamic 
Monte Carlo simulations employing embedded atom potentials for a 
model system. Three limiting models, all considering adatom diffusion, 
were employed to describe adatom deposition. The first model (A) con- 
siders adatom deposition on any free substrate site on the surface at the 
same rate. The second model (B) considers adatom deposition only on 
substrate sites which exhibit no neighboring sites occupied by adatoms. 
The third model (C) allows deposition at higher rates on sites present- 
ing neighboring sites occupied by adatoms. Under the proper condi- 
tions, the coverage(#) vs time(i) relationship for the three cases can 
be heuristically fitted to the functional form 9 = 1 — exp(—(3t a ), where 
a and f3 are parameters. We suggest that the value of the parameter 
a can be employed to distinguish experimentally between the three 
cases. While model A trivially delivers a = 1, models B and C are 
characterized by a < 1 and a > 1 respectively. 

Keywords: metal deposition, surface diffusion, dynamic Monte Carlo 
simulations 
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1 Introduction 



The study of the kinetics of formation and growth of two dimensional phases 
in electrochemical systems is a classical problem in the field of physical chem- 
istry. This sort of studies are experimentally performed by applying various 
perturbations to the system in order to obtain the response of a macroscopic 
observable which allows to infer the distinctive features of the nucleation 
and growth process. Over the last two decades, the improvement in in-situ 
nanoscopy techniques has also made possible to observe the time evolution of 
the surface morphology in electrodeposition experiments as it is done in the 
case of deposition experiments under ultra-high vacuum (UHV) conditions. 

From a theoretical point of view, a large variety of models are generally 
used to describe nucleation and crystal growth phenomena j2j. Many of these 
models apply continuum theories, while others use descriptions based on mi- 
croscopic concepts, as does the DDA (Diffusion, Deposition, Aggregation) 
model which will be referred to in section 4. Most of these models have been 
formulated to describe experimental situations similar to the ones set up un- 
der UHV conditions, where the adsorption process occurs randomly and at 
constant flux at any place on the surface. However, during electrodeposition 
of a metal, the entrance of ad-atoms to the system is a thermally-activated 
process in which the activation energy for ion reduction is determined mainly 
by the reorganization energy of the solvent This particular feature has 
raised a question about the type of site on the electrode surface where particle 
deposition occurs. It is possible that in some systems the reduction reaction 
occurs mainly on the edge of previously formed islands, whereas in other sys- 
tems particles deposit preferentially onto the clean terraces and then diffuse 
in order to join the growing islands as proposed by Bockris |3]. A conclusive 
answer to this question has not yet been found. This present work has two 
main aims. On one hand, we seek to study the most remarkable features 
of the electrochemical response and the evolution of surface morphology of 
the system for the two limiting cases mentioned above, in comparison with 
the trivial case where deposition occurs at the same rate on all free places of 
the surface. On the other hand, we shall give simple heuristic criteria that 
should allow the identification of the nature of the discharge process from 
experimental data. 

2 Approach to the problem 

As a model system, we consider here electrodeposition in the case of a neg- 
ligible adsorbate/substrate lattice misfit through the dynamic Monte Carlo 
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(DMC) method. In order to describe the interaction between particles in a 
realistic metallic system, we use semiempirical potentials corresponding to 
Ag deposition on Au(100). However, we believe that many features of the 
present model will also be valid for other hetheroepitaxy systems beyond the 
particular potential employed to describe the interaction between the par- 
ticles in the system. Furthermore, the DMC method allows to study the 
system evolution in real time. Concerning the type of perturbation applied 
to the system, in the present work we shall simulate a potentiostatic experi- 
ment, were the electrode potential is changed from an initial value where the 
surface of the electrode is free from adatoms, to a value where adsorption 
occurs. Furthermore, we shall consider irreversible deposition, that is, we 
neglect adatom dettachment from the surface. 

The relationship between the atom deposition rate on a given type of site 
x and the potential applied is given by the following equation[lJ: 

k dep>x (E) = k° dePtX a Me , + exp ^- A ^ p ' x ' j exp ^—2 — W 

where kd ep , x is the atom deposition frequency on site x , E is the measured- 
electrode potential vs. a reference electrode, AG^, is the activation energy 
for ion transfer from the solution to the crystal at E=0, a c is the charge trans- 
fer coefficient, om^+ is the activity of the metallic ion in the electrolyte, and 
k° depx is the rate constant for the deposition reaction on a site x. Accord- 
ing to equation (QJ), a simulation of adatom deposition appears a priori as a 
formidable task, since kd ep , x (E) should be calculated for every possible type 
of site on the surface. This would require a detailed knowledge of the elec- 
tron transfer itself. For example, the calculation of kd ep , x (E) could involve 
on one hand the local electronic properties of the site where the discharge 
of the ion is taking place, as well as the nature of the rearrangement of the 
solvent concomitant with the charge transfer process. However, as we stated 
in the introduction, in this work we shall be concerned with some special 
cases of equation (|T} and we shall think of the surface as essentially made of 
two types of surface sites: those which are next to at least one adatom and 
those which are not. Thus, in principle only two different values of kd ep , x (E) 
may occur. Among the various possibilities that emerge even after this sim- 
plification, the following three possible mechanisms for atom attachment to 
the growing phase were studied, which are considered in Figure 1: 

A Particle adsorption occurs on all unoccupied sites on the surface, which 
are considered as equivalent, independently of whether or not the site 
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is surrounded by atoms. In this case, a rate k is assigned to the pro- 
cess of atom entrance into an unoccupied site, independently of its 
surroundings. 

B Particle adsorption occurs only at sites corresponding to the terraces, 
i.e. where no adsorbed atoms are around. In this case an adsorption 
rate k is assigned to those sites. The adsorption rate is considered to be 
at those sites that have at least one of the nearest neighboring sites 
occupied, i.e. on the edge of steps or on kink sites. This model is based 
on Bockris' idea that the ion reduces preferentially on the terraces, 
since it loses there the least part of its solvation sphere. After that, the 
ion diffuses towards the step edges and then towards the kink sites. 

C Particles are considered to be discharged preferably at step edges. With 
this purpose, a rate ki is assigned to the entrance into terrace sites, i.e. 
where no neighboring atoms are around, and a rate &2 (where k% > k\) is 
given to the entrance of atoms into sites that have at least one atom as 
its nearest neighbor. This situation is comparable to the Avrami model 
[TTH |H] in which the monolayer grows from the island edges, assuming 
that atoms enter more easily there. 

In the present model the parameter k was changed between lCT 3 s -1 and 
10 +2 s _1 (at each adsorption site) and the nature of the site x was determined 
by its atomic enviroment as stated above. In the present model no assumption 
is made concerning the evolution of the surface morphology. This is a result 
of the deposition rate k and the interaction between the particles of the 
system that affects the diffusion properties. Thus emphasis is set on making 
a realistic model of atom diffusion on the surface, according to a potential 
that takes into account the many body effects and which is suitable for the 
study of metals. 

This paper is organized as follows: the model and the simulation tech- 
nique used are described in Section 3; the results are given in Section 4, 
which is subdivided into two parts: in the first, the system dynamic response 
is studied in terms of the evolution of the coverage degree and of the current 
(defined as i = d6/dt) as a function of time; in the second part, structural 
aspects are dealt with, that is, a description of the morphology and quantifi- 
cation of the formed islands is given. Finally, our conclusions are presented 
in Section 5. 
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3 Lattice Model and Simulation Technique 



3.1 Dynamic Monte Carlo Method 

Monte Carlo methods are used as computational tools in many areas of phys- 
ical chemistry. Although traditionally applied to obtain equilibrium proper- 
ties, they can also be used to study dynamic phenomena In order to do 
this, the following conditions must be fulfilled: a) in addition to satisfying 
the detailed balance criterion, the probabilities of transition must reflect a 
"dynamic hierarchy"; b) time increments between events must be correctly 
formulated in terms of the microscopic kinetics of the system; c) the events 
must be effectively independent. In the DMC method every step consists 
in a random selection of one of the possible processes. The probability of a 
process being selected is directly proportional to its rate. Once the randomly 
chosen method has been performed, all the possible processes are calculated 
and stored again in the corresponding vector and there is a time increment of 
At = —ln(u)/J2 v i where u is a random number between and 1, and J2 v i 
is the sum of the rates of all the possible processes. This time increment is 
due to the assumption that we are dealing with a Poisson process[5j. 

Silver adsorption on a defect-free Au(100) surface was studied. The com- 
putational model worked with a square arrangement of n x n adsorption sites, 
under periodic boundary conditions. The simulations were performed with 
n = 50 and n = 100. Unless otherwise stated, the results presented corre- 
spond to the former case. The sites corresponding to the nearest neighbors on 
each side were kept in a matrix. The system was initialized: for the present 
studies, the initial state corresponded to a clean Au(100) surface, but the 
model could be easily extended to surfaces exhibiting previously adsorbed 
silver or gold atoms. The rates of every possible process were stored in a 
vector to perform the random selection described above. In this work since 
we consider the case of irreversible deposition (no desorption), such processes 
were the entrance of an atom into an empty site and the motion of an atom 
from a site to one of the four sites corresponding to first nearest neighbors 
(face 100). For example a fixed rate may be assigned to the first process, in- 
dependently of the occupation of the neighboring sites. This would cause an 
increase in the coverage degree given by a law of the type 9=1 — exp(-kt), 
where t denotes the time elapsed. In the most general case, the entrance 
rate of a particle depends on whether or not the site is surrounded by other 
atoms, as discussed below. 

The algorithm developed by Hoshen and Kopelman was used to calculate 
the island number and size [9] and average was taken over 10 simulations for 
each rate. The interval was divided into 100 bins in order to calculate the 
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average at each coverage degree. 

3.2 Energy Calculation 

To calculate the activation energies for adatom diffusion the Embedded Atom 
Method (EAM) was used This method takes into account many body 
effects; therefore, it represents better the metallic bonding than a pair po- 
tential does. The total energy of the system is calculated as the sum of the 
energies of the individual particles. Each energy is in turn the sum of an 
embedding (attractive) energy and a repulsive contribution that arises from 
the interaction between nuclei. The EAM contains parameters which were 
fitted to reproduce experimental data such as elastic constants, enthalpies of 
binary alloys dissolution, lattice constants, and sublimation heats. We have 
employed the EAM within a lattice model as described in references 0IH|- 

3.3 Diffusion Rate Calculations 

Diffusion rates used in the simulation were previously tabulated according 
to the different environments that an atom can find on the surface. Various 
configurations were generated, consisting of different arrangements of atoms 
near the starting and ending sites. For each configuration the path followed 
by an atom to jump from a site to the neighboring one was traced and 
the energy in each position was calculated, minimizing it with respect to 
z (coordinate perpendicular to the surface plane). The activation energy 
E a was calculated as the difference between the saddle point and the initial 
minimum in the energy curve along the reaction coordinate. The vibrational 
frequency v of the atom in the starting site was calculated performing the 
harmonic approximation near the minimum of the curve. The diffusion rate 
was calculated as v = v exp(—E a /kT). More details of this calculation 
procedure can be found in reference |Bj 

4 Results 

4.1 Evolution of the coverage degree and the current as 
a function of time 

4.1.1 Model A: Adsorption Independent of Surroundings. Dy- 
namic response. 

Model A is the simplest one and the evolution of the coverage degree responds 
to the equation d9/dt=k(l-9), where k denotes the adsorption rate per site. 
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In this case 9=l — exp(—kt) and the current satisfies the law i = k exp(-kt). 
In a plot of the coverage degree as a function of kt, the curves corresponding 
to the different simulations overlap, i.e. the evolution of 9 as a function of kt 
is independent of the atom entrance rate. However, as we shall later see, the 
surface structure differs significantly regarding the island number and size. 
Note the difference of the present electrochemical (potentiostatic) conditions 
with respect to the UHV ones. In the latter case, the coverage degree evolves 
linearly with time, while in the potentiostatic case the rate is proportional 
to the available surface resulting in an exponential law. 

4.1.2 Model B: Preferential Adsorption on Terraces. Dynamic 
response. 

For case B, that is, with atom adsorption only on the terraces, Figure 2 shows 
the evolution of 9 as a function of kt for several k values. Here it can be 
observed that for higher entrance rates the surface tends to a lower coverage, 
i.e. it takes longer to reach a high coverage, while for lower rates the coverage 
tends to be closer to one. This is so because at higher rates more and smaller 
islands are formed, generating a larger amount of edge sites where deposition 
is not allowed. On the other hand, at low deposition rates the particles have 
enough time to diffuse reaching a more compact distribution, i.e. forming 
fewer and larger islands, therefore leaving more sites available for adsorption. 
Further analysis shows that at relatively short times 9 initially evolves like 
9 = 1 — exp(—j3t a ), with a<l. This can be verified by a plot where ln(-ln(l- 
9)) is represented as a function of ln(t) (Figure 3), since straight lines of slope 
a and y-intercept lnf3 are obtained. Table 1 shows a and (5 values for the 
various k. For higher rates we found a values considerably lower than those in 
case A (equivalent sites) where a= 1. For lower rates a increases tending to 
1. If the system is allowed a longer time to evolve, a change in the exponent 
is observed, i.e. the curve becomes much more horizontal. This effect is much 
more pronounced at high adsorption rates. The physical explanation for this 
is that at a given moment many small islands have already formed, and 
consequently few unoccupied sites are left on the terrace, which makes the 
entrance of new atoms more difficult. At this point diffusion starts playing a 
fundamental role, since the islands start rearranging in order to form larger 
ones. This process is slow compared with the diffusion of free atoms on the 
terrace. As diffusion occurs, adsorption sites are generated, allowing a slow 
entrance of new atoms, which produces a small current. Finally, at high 
coverage vacancy islands and simple hollows are left . There the diffusion 
rate of vacancies is higher and so is the generation of adsorption sites, which 
causes the curve slope to increase at the final stage. 
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Figure 4 shows the evolution of 6 and of the current as a function of time 
for k = 1 s -1 . Here a comparison can be made between model A (adsorption 
equal at all sites), and model B (adsorption only on the terraces). In the case 
of the latter, the coverage degree is found to be lower and tends to 1 more 
slowly for the reasons stated before. For the former case, current decreases 
exponentially, whereas for adsorption only on the terraces, it decreases faster 
at the beginning and more slowly afterwards, satisfying the law dO/dt = 
fia.t a ~ x exp(— (3t a ) . The points obtained as the result of taking average over 
10 simulations are shown as well as the curves predicted by this ansatz, 
generated with the values of a and (3 obtained from the fit. 

4.1.3 Model C: Preferential adsorption on the Edges. Dynamic 
response. 

The third case corresponds to atom adsorption that is preferential on step 
edges. Simulations with different adsorption rates on the terraces (ki) and 
on the edges (k 2 ) were performed. Table 2 shows the rate constants of all 
the systems that have been studied here. 

In these simulations the time at which the first particle deposits is taken 
as the initial time, since the time elapsed until this first process occurs may 
be, in some cases, significantly longer than the time taken by the simulation, 
especially when particle adsorption occurs on the edges of islands at rates 
much higher than those on the terraces. Figure 5 shows the evolution of 
the coverage degree as a function of time for the cases having k 2 = 1.0s" 1 
as entrance rate on the step edges, and k\ = 0.01s -1 , ki = 0.001s" 1 and 
ki = O.OOOls^ 1 as entrance rates on the terrace. In these three cases, as well 
as in all cases where the adsorption rate is higher on step edges than on the 
terrace, it can be noticed that the initial evolution of the coverage degree has 
a curvature that first points upwards and then downwards, i.e. it presents 
a saddle point. This is due to the fact that the current depends mainly on 
the amount of edge sites, which initially increases and afterwards decreases 
as coalescence occurs. 

In these before, the coverage degree was assumed to satisfy a 

law 9 = 1 — exp(—f3t a ) and was represented in a plot as ln(— ln{l — 6)) as a 
function of ln(t — t ) in order to obtain the coefficients a and (3. These are 
summarized in Table 3. These values give a qualitative idea of the tendencies, 
since they correspond to a single simulation. It can be observed that, in all 
cases, a is larger than 1, and it tends to decrease for a given value of k\ as ki 
decreases. Figure 6 shows the current as a function of time, along with the 
curves d6/dt = /3cd a ~ 1 exp(— (3t a ), with the parameters obtained from the fit. 
Three examples are shown. The plot at the bottom of the page corresponds 
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to k\ = 10 _4 s _1 and k 2 = 10 2 s _1 . Here the two rates are very different from 
each other and it can be clearly observed that initially the current is 0, then 
it increases up to a peak , and after that decreases. The plot in the middle 
corresponds to an intermediate case, where the rates are k\ = 10 _3 s _1 and 
k 2 = 1.0s" 1 and the curve is qualitatively very similar to that of the previous 
case. The plot at the top corresponds to k\ = lCT 2 s _1 and k 2 = 2 x 10 _2 s _1 . 
In this case the two rates are very similar and a current maximum is not 
observed in the simulated curves. This case is more similar qualitatively to 
the exponential decrease corresponding to equal adsorption at all sites. 

4.2 Structural Aspects 

In this section we shall mention first some of the quantities often used to 
describe the dynamic behavior of the surface structure in the so-called atom- 
istic crystal growth models. We shall especially make use of the concepts 
used in the DDA model, which considers three basic processes as responsible 
for epitaxial growth: Deposition, Diffusion, and Aggregation [HI [T2j . Then, 
we shall study the behavior of these quantities in each of the three cases con- 
sidered in this work and their relationship with the corresponding current 
transients described in the previous section. 

4.2.1 Use of the Island-Size Distribution Function to Describe 
Surface Morphology 

A fundamental quantity in describing kinetically the growth of epitaxial 
monolayers and submonolayers is the island-size distribution function N s (t), 
which is the density per site of islands of size s at time t, s being the number 
of atoms in the island. The function N s (t) is used in the DDA model to 
describe the surface structure and it has mainly been applied to the analy- 
sis of experimental results obtained under UHV conditions. In this case the 
coverage degree 9 changes linearly with time and then the function of size 
distribution is generally expressed as a function of 6. In this work we shall 
follow the same criterion, though we should keep in mind that in our simula- 
tions 6 does not evolve trivially, because of which many of the relationships 
found do not strictly correspond to those of the dynamic scaling. 

Defining the total number of islands N and the coverage degree 6 for a 
given instant by 

N = Y,N S 6 = J2sN s (2) 

s>2 s>l 
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then the average island size S(8) can be written as a function of the zeroth 
and first moments (J2 s >2^s, J2 s >2 s ^s) as 

m - y^K - — (3) 

In analyzing the evolution of the surface structure and its relationship to 
the evolution of the current and of the coverage degree, we shall distinguish 
various dynamic regimes in the behavior of the total island density N and of 
the monomer density Ni. Then, we shall present our work in a way similar 
to the one by Family et al [TT] . 

It is important to highlight that the scaling theory is satisfied only under 
certain deposition conditions (e.g. relatively low deposition rates), and nat- 
urally, only during the aggregation regime. Thus only some of the performed 
simulations can be expected to satisfy completely the observations predicted 
by this theory. Such is the case for models A and B when the coverage degree 
is approximately at the interval (0.1,0.4) and the deposition rate is low (low 
overpotentials) , as we shall see in sections 4.2.2 and 4.2.3. 



4.2.2 Model A: Adsorption Independent of Surroundings. Struc- 
tural Aspects 

As mentioned in Section 4.1.1, in model A the evolution of the coverage 
degree and of the current are given trivially by an exponential function. 
However, the observed variety of surface morphology depending on the de- 
position conditions indicates that no relationship exists between the surface 
structure and the electrochemical response. 

Six different entrance rates of atoms were studied: 10 +2 s -1 , 10 +1 s -1 , 
10°s~ l , 10 _1 s~ 1 , 10~ 2 s _1 and 10 _3 s _1 . Figure 7 shows the evolution of the 
surface for two different adsorption rates: the left-hand column corresponds 
to k — 10 +2 s _1 , and the right-hand column to k = 10~ 2 s -1 . From top to 
bottom the coverage degrees are 6 = 0.1, 9 = 0.3 and 9 = 0.6. As it can 
be observed, at higher adsorption rates more and smaller islands are formed. 
This is true when coverage is low or intermediate, since at high coverage and 
rates islands are branched so they coalesce in order to form one island. At 
high rates and low coverages a large number of monomers can be observed. 
This is due to the fact that, when particles have such a high entrance rate, 
they have little time to diffuse, whereas at a low entrance rate particles have 
time to join the already existing islands. 

Figure 8 shows the total island number, and the monomer number divided 
into the total site number for the various entrance rates of atoms. The 
upper and the lower part of the figure corresponds to the island density 
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and to the monomer desity respectively. As it can be observed, at higher 
deposition rates both the monomer density and the island density are larger. 
In all cases the monomer maximum is reached at lower coverages than the 
maximum in the island density. At low deposition rates three zones can 
clearly be distinguished during deposition. The first one consists in an initial 
increase in the island number, which corresponds to the nucleation process. 
Nuclei are formed at the beginning, so we can refer to an instantaneous 
nucleation. Then, the number of islands remains essentially constant and 
this corresponds to the aggregation regime in which every adatom joins an 
existing island whose size increases in this way, but whose number does not. 
Finally, the coalescence regime sets in, during which the number of islands 
decreases as they begin joining one another to form larger islands. Eventually, 
only one island is formed which grows up to form the complete monolayer. 
At higher deposition rates, the island density reaches a maximum at higher 
coverage degrees, and we can refer to a progressive nucleation, since the 
number of islands is increasing over a larger zone of 9 and it is not possible 
to differentiate a specific aggregation stage. 

Figure 9 shows the average size of islands (calculated as the average num- 
ber of atoms of every island divided into the total number of sites) as a 
function of the coverage degree for the various adsorption rates. In all cases 
there is, initially, a linear increase in the island size, with higher slopes for 
lower rates. At higher rates, there is a very small increase in the average 
island size, since new small islands are continuously forming. On the other 
hand, at low rates few islands are formed, whose size increases considerably 
as new atoms join them. After that there is an abrupt change in size, which 
corresponds to coalescence, and the last stage observed corresponds to size 
increase of a single island. 

Percolation properties of the system were also analyzed. The percolation 
coverage degree 9 P is defined as the coverage degree at which at least one 
island is formed that crosses over the system, at least in one direction. Figure 
10 shows 6p as a function of log(R), with R — D/k, where D is the diffusion 
rate of the atoms on the surface (D = 436 in this case). This was 
calculated for a 100 x 100 surface and averaged over several simulations. 
As it can be observed, 9 P increases as R does, that is, as atom adsorption 
rate decreases. This is due to the fact that at lower rates islands are more 
compact, therefore they take more time to join one another. On the other 
hand, at high adsorption rates islands branch and they are more likely to 
join together to form a single island at lower 6*'s. 
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4.2.3 Model B: Preferential Adsorption on Terraces. Structural 
Aspects. 

As it was previously detailed, model B considers that adatom deposition 
occurs only at sites with coordination 0. In electrochemical terms, this means 
that the reduction reaction occurs only on defect-free terraces. In this case, 
we have shown above that the evolution of the coverage degree (and of the 
current) is non-trivial and is given by a "stretched exponential" function 
which emerges from a scenario similar to that proposed by Palmer et. al. 
[ISlIIlj. The exponent a depends on the applied overpotential, and in all 
cases is less than one. Simulations performed under conditions analogous to 
those of Figure 7 yielded images that are qualitatively very similar to those 
of model A presented there. However, it was found that for model B the time 
at which each degree of coverage occurs is longer, that is to say, the filling in 
of the surface is slower. Comparative figures are given in Table 4. It can be 
noticed that this effect is stronger at higher coverages and at high deposition 
rates. As explained before this is due to the lesser availability of adsorption 
sites. 

In Figure 11 we can observe the evolution of the total density of islands 
N and the density of monomers N\ as a function of 6. A comparison between 
this plot and the one of the previous model in Figure 8 shows that at high 
rates the island and monomer density are larger, while at low rates the curves 
practically agree. This is due to the fact that at high adsorption rates atoms 
have little time to diffuse and reorder, and in turn those that adsorb do 
it in sites not surrounded by atoms. All this causes the formation of a 
larger amount of islands, whereas at low rates, in spite of adsorption only at 
neighbor- free sites atoms have sufficient time to diffuse and join the existing 
islands and the net growth occurs in a way similar to that of the previous 
model. In addition to the larger island density we can observe a shift at high 
coverage (between 0.4 and 0.6), which indicates that there are more islands 
in that zone, that is to say, the islands take longer to coalesce and do it at 
higher 9s. 

Figure 12 shows the average size of islands as a function of the coverage 
degree. If we compare this plot with the one of Model A we can observe 
that the abrupt change in size shifts to higher 9s, especially for the case of 
an atom entrance rate of 10 _3 s _1 , which is the slowest rate. This is due 
to island compactness, which causes islands to take longer in joining and 
forming a single one. 
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4.2.4 Model C: Preferential Adsorption on the Edges. Structural 
Aspects. 

This model considers that atoms on island edges adsorb at a rate higher than 
that of atoms which adsorb on the terrace. Adsorption rates of 10 _2 s _1 , 
10~ 3 s _1 and 10~ 4 s _1 on the terrace (ki), and up to 10 +2 s _1 on the edges(/c 2 ) 
were studied. For very large rate differences a single island is formed in our 
simulation system, which grows from the edges. As adsorption rates become 
similar the formation of more islands is observed. Figure 13 shows the state 
of the surface at 9 = 0.1 for different values of k\ and k 2 , entrance rate on 
the terrace and on the edges respectively. The left-hand column corresponds 
to k 2 = 10 +2 s _1 . Under these conditions, only for k\ = 10~ 2 s~ 1 two islands 
were formed, while at lower values of k\ a single island appears. On the other 
hand, the results on the right column show that for similar values of k\ and 
&2 several islands are formed. In all cases, however, the island shape is quite 
compact and dendritic shapes are not observed, as they were found at high 
rates in the other models. 

5 Conclusions 

In this paper we considered several aspects of 2-D metal deposition by means 
of dynamic Monte Carlo simulations As a model system we considered Ag 
deposition on Au(100) under various deposition conditions using interatomic 
potentials adequate for metals and a realistic diffusion simulation model. 

The reduction of Ag + ions was modeled in three different ways: in the 
first and simplest case, deposition of atoms was assumed to occur on any 
unoccupied site of the lattice independently of the site surroundings. Here 
the evolution of the coverage as a function of time was given in a trivial way 
by an exponential function, while the structure of the surface depended on 
the deposition rate. In the second case, particle adsorption was allowed only 
on defect-free terraces, that is, only on those sites characterized by adsorbate- 
free surroundings. It was observed that the evolution of the coverage rate is 
given by a function 9=1 — exp(— f3t a ) , where the exponent a depends on the 
deposition rate, that is, on the overpotential applied, though in all the cases 
considered it turned out to be smaller than one. At high deposition rates 
a change in this exponent was also observed, which indicates different types 
of evolution of the surface structure. Thus a relationship exists between the 
coverage evolution - and therefore the current evolution - and the change in 
the surface structure. In the third deposition model, the ion reduction was 
considered to be possible on edge sites as well as on sites located on clean 
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terraces. A higher deposition rate was always used on the edges. In this 
case, the coverage degree seems to follow again a 6 = 1 — exp(—j3t a ) law, 
but with exponent a always larger than one. This kind of behavior can be 
compared with the predictions of the Avrami model |3], where the exponent 
a is two or three, depending on whether the nucleation is progressive or 
instantaneous. As regards the evolution of surface morphology, there was a 
remarkable reduction in the number of islands, and a single island was formed 
in our simulation system when there was a very large difference between 
the rates over the two types of sites. The exponent values reported must 
be considered to give a qualitative tendency, in the sense that we did not 
study systematically the finite-size effects and they correspond to a single 
simulation run. These features may be important in determining the exact 
values of the exponents. 
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7 Tables 



Table 1 

Fitted a and f3 values from simulations with various adsorption rates k 
onto an unoccupied site . These parameters were obtained by fitting of the 
law 9 = 1 — exp(—j3t a ) from simulations with model B described in the text. 



k(rate) 


a 


P 


10 2 


0.722 


12.72 


10 1 


0.748 


2.914 


10° 


0.799 


0.588 


KT 1 


0.868 


0.095 


io- 2 


0.909 


0.012 


1(T 3 


0.951 


0.001 



Table 2 

Parameters employed in the simulation studies with model C. k\ and &2 
denote the adsorption rates on sites without adsorbed neighbors( 'terrace site 
') and with adsorbed neighbors( 'border site ') respectively. Crosses indicate 
the sets of values of rate constants considered . 



k 1 (terrace)\k2 (border) 


10 2 


10 1 


10° 


10" 1 


2 x 10~ 2 


io- 2 


2 x 10~ 3 


io- 2 


X 


X 


X 


X 


X 






10~ 3 


X 


X 


X 


X 




X 


X 


10" 4 


X 


X 


X 


X 




X 





Table 3 

Fitted a and (3 values for various simulation conditions with model C. 
They are obtained by fitting of the law 9 = 1 — exp(—j3t a ) . 
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kl(terrace) 


k2(border) 


a 





io- 2 


IO 2 


2.3 


63 


10- 2 


10 1 


2.4 


1.3 




10° 


1.9 


0.03 


IO- 2 


io- 1 


1.5 


0.005 


IO- 2 


2 x IO" 2 


1.2 


0.006 


io- 3 


10 2 


2.1 


36. 


10- 3 


10 1 


2.0 


0.51 


IO" 3 


10° 


2.0 


0.009 


IO" 3 


io- 1 


1.7 


7.1 x 10- 4 


IO" 3 


IO" 2 


1.3 


4.3 x 10~ 4 


IO" 3 


2 x IO" 3 


1.1 


7.4 x 10" 4 


10~ 4 


IO 2 


2.0 


27. 


10~ 4 


10 1 


1.9 


0.24 


IO' 4 


10° 


2.0 


0.002 


io- 4 


io- 1 


2.1 


2.9 x 10- 5 


10~ 4 


IO" 2 


1.7 


4.9 x IO" 6 
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Table 4 Time t$ at which a given degree of coverage 9 occurs in models 
A and B described in the text, k denotes the rate assigned to the process of 
atom entrance onto an unoccupied site. 



e 


t e (k = 10 2 s- r ) 


t (k = 10 2 s- r ) 


t e (k = 10-' 2 s- r ) 


t e (k = W^s- 1 ) 




mod el A 


model B 


mod el A 


model B 


0.1 


0.00094 


0.0013 


11.3 


11.6 


0.5 


0.007 


0.020 


65.5 


88.2 


0.7 


0.012 


0.083 


114. 


177. 


0.9 


0.024 


194. 


221. 


430. 



8 Figure Captions 

Figure 1: Scheme of the mechanisms for atom attachment to the growing 
phase 2-D studied in the present work: 

A) Particle adsorption occurs on all unoccupied sites on the surface at the 
same rate k. 

B) Particle adsorption occurs only at sites where no adsorbed atoms are 
around ( 'terrace sites ') at the rate k. The adsorption rate on sites 
that have one or more neighboring adatoms is 0. 

C) Particles are considered to be discharged preferably at sites with neigh- 
boring adatoms ( 'edges '). A rate k\ is assigned to the entrance onto 
terrace sites and a rate k 2 ( k 2 > k±) is given to atoms entrance into 
sites that have at least one neighboring adatom. 

Figure 2: Evolution of the coverage degree 9 as a function of k * t corre- 
sponding to the adsorption model B described in the text ('terrace ' adsorp- 
tion). The adsorption rate k is given in the figure in units of s -1 . 

Figure 3: Plot of ln{— ln(l — 9)) as a function of ln(t) at various deposi- 
tion rates k for Model B described in the text ('terrace ' adsorption). The 
adsorption rate k is given in the figure in units of s -1 . Slopes (exponents) 
associated to the linear relation found during the first stage of deposition are 
detailed in table 1. 

Figure 4: Evolution of the coverage degree 9 and of the current for models 
A and B for a deposition rate of 1 s _1 . 
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Black solid line: evolution of 6 in model A. Gray dashed line: evolution 
of 9 in model B. 

Circles: current evolution in model A. The gray solid line represents the 
theoretical fitting according to the equation d6/dt = /3at a ~ 1 exp(—j3t a ) 

Diamonds: current evolution in model B. The black dashed line represents 
the theoretical fitting according to the equation d6/dt = fioct (X ~ 1 ex'p{— /3t a ) 

The exponents a and (3 are given in Table 1. 

Figure 5: Evolution of 9 as a function of time for the deposition model C. 
In these simulations, a deposition rate k 2 = 1 s _1 is considered on the edges 
and different deposition rates k\ are assumed at the terraces. The k\ values 
are given in the figure in units of s~ x . 

Figure 6: Current (d9/dt) as a function of time for deposition model C. 
The full line are simulation results, and the dashed lines show the curves 
predicted by the function d6/dt = j3at a ~ l exp(— f3t a ) . Each plot employs the 
parameters obtained through minimum square fits of the simulation data. 
The adsorption rates employed for the simulations were the following: Top 
box: h = lO^s- 1 ;k 2 = 2x lCT 2 ^ 1 ; middle box: h = lCT 3 ^ 1 ; k 2 = Is -1 ; 
box at the bottom: k\ = 10~ 4 s _1 ; k 2 = 10 2 s _1 . 

Figure 7: Frames showing surface morphology in model A at three dif- 
ferent values of 6. The deposition rates were k = 10 2 s _1 for the left-hand 
column and k = 10^ 2 s~ 1 for the right-hand column. The coverage degrees 
are, from top to bottom, 9 = 0.1, 9 = 0.3 and 9 = 0.6. 

The associated times are t = 9.4 x 10~ 4 s, t = 3.5 x 10~ 3 s and t = 
9.2 x 10 _3 s from top to bottom for the left column and t = 11.31s, t = 34.22s 
and t = 87.26s for the right one. 

Figure 8: Plot of island density (top) and monomer density (bottom) as 
a function of the coverage degree for different adsorption rates according to 
model A. Rates used in each case are indicated in the body of the figure in 
units of s -1 . 

Figure 9: Average island size < s > (defined as the average number of 
atoms per island divided by the total number of sites iV^ ) as a function of 
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the coverage degree 9, for different adsorption rates according to model A. 
The rates used in the simulations are included in the figure in units of s^ 1 . 
N T = 2500. 

Figure 10: Coverage degree at which percolation occurs, 9 P as a function 
of log(R), where R = f , D being the diffusion rate of an adatom on the 
clean surface and k the particle entrance rate. Simulations were performed 
on a 100 x 100 square lattice. 

Figure 11: Plots of island density (top) and monomer density (bottom) 
as a function of coverage degree for different adsorption rates according to 
model B. Rates used in each case are indicated in the body of the figure in 
units of s _1 . 

Figure 12: Average island size < s > as a function of the coverage degree 
9, for different adsorption rates according to model B. The rates used in the 
simulations are included in the figure in units of s _1 . 

Figure 13: Frames showing the surface structure in simulations with 
model C at a coverage degree 9 = 0.1 for different adsorption rates on ter- 
races and edges. Frame 1: k± — 10 _2 s _1 ; k 2 = 10 2 s _1 ; Frame 2: k\ = 
Kr 2 s _1 ; k 2 = 2x 10- 2 s _1 ; Frame 3: h = lO^s' 1 ; k 2 = 10 2 s~ 1 ; Frame 4: 
ki = lO- 3 ^- 1 ; k 2 = 10- 2 s~ 1 ; Frame 5: h = KTS -1 ; k 2 = 10 2 s" 1 ; Frame 6: 
k ± = lO^s" 1 ; k 2 = 10- 2 s-\ 
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